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For Int. J. Quantum Chem. 

This note summarizes the motivation for extending current 
density-functional theory to include nonlocal one-electron po- 
tentials, and proposes methodology for practical calculations. 
The theoretical model, orbital functional theory, has been 
shown to be exact in principle for the general N-electron prob- 
lem, but must depend in practice on a parametrized correla- 
tion energy functional. The discussion here is intended to 
honor Lee Allen and to bring up to date some conversations 
that we began in 1954. 



I. INTRODUCTION 

The density functional theory (DFT) of Hohenberg, 
Kohn, and Sham [1,2] has become the method of choice 
for computation of electronic structure and properties 
of large molecules and of solids. This theory justifies 
an independent-electron model that is considered to be 
exact in principle for N-electron ground states [3,4]. It 
is widely assumed that such an exact model can be ex- 
pressed in terms of local potential functions for exchange 
and correlation, although evidence to the contrary has 
existed for some time [5]. It has only recently been 
shown that this locality hypothesis is inconsistent with 
exact theory for more than two electrons [6] . The specific 
mathematical point is that density functional derivatives 
of the Hohenberg-Kohn universal functional cannot be 
equivalent to local potential functions, as assumed in es- 
tablished literature [3,4], for more than two electrons in 
a compact electronic system. These functional deriva- 
tives can be constructed explicitly [6,7], and are found to 
be Gateaux derivatives, the generalization to functional 
analysis of analytic partial derivatives [8] . For more than 
two electrons, they do not reduce as commonly assumed 
to simpler Frechet derivatives, equivalent to multiplica- 
tive local potential functions, the generalization to func- 
tional analysis of analytic total derivatives [8] . 

This mathematical development explains a number of 
inconsistencies or paradoxes that result from assuming 
the locality hypothesis [5] . Failure of locality for the ex- 
change potential in a Hartree-Fock variational model of 
DFT is implied by results of optimized effective potential 
(OEP) calculations going back as far as Aashamar, Luke 
and Talman [9], confirmed by recent quantitative tests 
[10]. The resulting problems or paradoxes include fail- 
ure of DFT linear response theory [11] to agree in the 
exchange-only limit with the time-dependent Hartree- 
Fock theory of Dirac [12]. This failure can be traced 
to the locality hypothesis [13]. 



The most striking paradox affects the kinetic energy. 
In contrast to the locality hypothesis, Kohn and Sham 
use the kinetic energy operator of Schrodinger for what 
should be the density functional derivative of a ground- 
state kinetic energy functional. If the locality hypothe- 
sis were a rigorous consequence of variational theory, it 
should be valid for the kinetic energy functional, imply- 
ing for noninteracting electrons an exact Thomas-Fermi 
theory (TFT) equivalent to DFT. However, the varia- 
tional equations of these theories are incompatible for 
more than two electrons [14]. The Kohn-Sham equations 
of DFT are correct for atomic shell structure, but the lo- 
cality hypothesis in TFT is inconsistent with the exclu- 
sion principle [14]. This is evident for the lowest ls2s 3 S 
state of an atom with two noninteracting electrons, per- 
haps the simplest unequivocal example of this inconsis- 
tency. Unless both Is and 2s orbital densities are inde- 
pendently normalized, the lowest total energy is achieved 
for two electrons in the Is state. Fermi-Dirac statistics 
require an independent Lagrange multiplier ei for each 
orbital, while TFT provides only a single parameter fj, 
for the total density. An equivalent inconsistency is im- 
plied for interacting electrons. This conclusion has been 
disputed and defended in recent publications [15,16,7]. 

Given a choice between an exact theory and one that 
can be expressed entirely in terms of local potential func- 
tions, there are pragmatic arguments supporting either 
alternative. The existing computational methodology of 
DFT requires use of local potential functions. In the 
local-density approximation (LDA), which is entirely jus- 
tified by variational theory [7], first-principles calcula- 
tions for condensed matter have had a revolutionary im- 
pact on understanding physical properties due to elec- 
tronic structure. However, for magnetically ordered ma- 
terials, there is a persistent problem due to the non- 
physical electronic self-interaction, which may require use 
of nonlocal potentials in a fully predictive theory. For 
molecules, long-range potentials due to polarization re- 
sponse can become significant in determining molecular 
conformations, but apparently cannot be modeled using 
only approximate local correlation potentials. The opti- 
cal potentials required in electron scattering theory are 
not directly described by energy-independent local po- 
tentials. 

With such applications or extensions of existing theory 
in mind, the present paper will consider practical issues 
that arise in extending existing local-potential method- 
ology to the nonlocal potential functions appropriate to 
an ultimately exact independent-electron model. 
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II. VARIATIONAL THEORY 

The variational theory of independent-electron models 
is most simply developed as an orbital functional the- 
ory (OFT) [17]. For an N-clcctron cigcnstate such that 
(H - E)^ = and any rule * — > $ that determines a 
model or reference state <£>, unsymmetric normalization 
($|*) = ($|$) = 1 implies that E = = E +E c , 

where Eq = ($|iJ|<!>) is an explicit orbital functional, 
and E c = (<I>|.ff"|\f' — $) defines the correlation energy. 
Restricting the discussion to nondegenerate states, $ 
is a Slater determinant constructed from occupied or- 
bital functions {<pi} 1 with occupation numbers rn = 1. 
Spin indices and sums are assumed here but are sup- 
pressed in the notation. For orthonormal orbital func- 
tions, E = T + U + V, where 



5F 



T = Y,m(i\i\i); V = ^2m(i\v\i), 



(1) 



for t = —\y 2 - Two-electron functional are defined for 
u = l/r 12 by U = E h + E x , where 

Eh = o X! njn ^H^) ; Ex = ~o ^ninjivHji)- ( 2 ) 



If Q = I - $$ f , E c = - $) = ($|i?|Q*). This 

implies an exact but implicit orbital functional [17] 



E c = -(<b\H[Q(H -E -E c - i^Q^H^), 



(3) 



for r\ — ► 0+. 

Orbital Euler-Lagrange (OEL) equations follow imme- 
diately from standard variational theory [18], in terms of 
the orbital functional derivatives 
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(4) 



using u = v h (r) + v x , where is the classical Coulomb 
potential, and v x is the Fock exchange operator. For 
independent subshcll normalization, Lagrange terms 
[{i\i) — l]ei are subtracted from the energy functional. 



Defining Q(pi 
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, the variational equation is 



J d 3 rJ2ni(S<t>*{G-£i}<f>i + cc) = 0, 



(5) 



for unconstrained orbital variations in the usual Hilbert 
space. This implies the OEL equations 



{Q - ei}<f>i = 0, i= N. 



(6) 



Orbitals of different energy are orthogonal. Equivalently, 
the universal functional F = E — V defines 



- = T4>i = {i +u + v c }(f>i. 



and the OEL equations are 



Tfyi = {e 8 - v}<f)i, i = 1, • • • , N. 



(7) 



(8) 



Overlap products in ($|i7|$) drop out of this argu- 
ment, because they either produce a unit factor or arc 
multiplied by a factor that vanishes. Explicitly, for non- 
interacting electrons, the same OEL equations are ob- 
tained if the simplified functional Eo = T + V is replaced 
by ($|f + V\$). The functional differential 

J^riif d 3 r{Scf>*(r){t + v(r) - e 4 }^(r) + cc]II j¥i (j|j) 

+ J2n i (i\t + v(r)-e i \i)SU j ^ i (j\j). (9) 

i 

must vanish for unconstrained orbital variations about a 
stationary state. It cannot do so unless all Lagrange mul- 
tipliers ti are included, precluding an equivalent Thomas- 
Fermi theory. Variations of the overlap product are mul- 
tiplied by a factor that vanishes. 

III. THE NEED FOR NONLOCAL POTENTIALS 

Hohcnbcrg-Kohn theory can be proven for a nonde- 
generate ground state of any version of OFT in which 
E xc = E x + E c is approximated by a parametrized func- 
tional of p [19]. Model electronic density p = J2i n iPi 
is a sum of orbital subshell densities Pi(r) = 4>*(r)<fti(r). 
Defining v xc (j)i — SE xc /md(j)* , the corresponding Kohn- 
Sham theory is determined by the OEL equations given 
above. The still controversial issue of whether or not the 
effective one-electron potentials in these equations can be 
expressed as multiplicative local potential functions de- 
pends on the relationship between the orbital functional 
derivatives of Eqs.(4) and the density functional deriva- 
tives needed for DFT. 

Density functional derivatives must be consistent with 
the OEL equations. A crucial logical point is that be- 
cause functional definitions must be extended to include 
variations that are not constrained by normalization 
[18], variational theory restricted to normalized ground 
states [20] cannot determine OEL equations. A varia- 
tional derivation of the Kohn-Sham equations must ex- 
tend the universal ground-state functional of Hohenberg 
and Kohn to density variations driven by unrestricted 
variations in the orbital Hilbert space. This extended 
functional is denoted here by F s . Orbital theory simpli- 
fies the mathematics, because such variations are neces- 
sary and sufficient to determine the Euler-Lagrange equa- 
tions [18]. 

The rigorous theory of Englisch and Englisch [20] 
proves the existence of a generic density functional 
derivative for a nondegenerate ground state. Although 
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not considered in the original derivation, it is clear from 
the OEL Eqs.(8) that the undetermined constant in 
[20], Eq.(4.1), must have orbital indices in order to be 
consistent with the electronic subshell structure implied 
by Fermi-Dirac statistics [7]. This is made explicit by 
Eqs.(7, 8) here, which imply that a functional derivative 
exists that depends on an orbital subshell index. Evalu- 
ated for a ground state [6,21], the functional differential 

SF S = / rf3r iM ( r )^"<M r ) + cc} 



= Yl Ui J rf M e i ~ v(r)}5pi(r) 



determines the functional derivative 



SF S 
riiSpi 



v(r). 



(10) 



(11) 



Because this depends on an orbital index, it defines a 
Gateaux derivative [6,8]. If a Frechet derivative were to 
exist, equivalent to a multiplicative local potential func- 
tion, it would imply the Thomas-Fermi equation, 



SFs^ 
Sp 



= p - v(r). 



(12) 



Because all dF s /riiSpi must equal SF s /5p if the latter 
exists, these equations are inconsistent unless all m are 
equal. Hence a Frechet derivative cannot exist for more 
than one electron of each spin without violating the ex- 
clusion principle. This result is not in conflict with rig- 
orous analysis [20], which establishes the existence of 
Gateaux derivatives in general, but cannot distinguish 
between Gateaux and Frechet derivatives because nor- 
malization is constrained [6]. For nonintcracting elec- 
trons, this analysis applies to the Kohn-Sham kinetic en- 
ergy functional T s [7] . If independent subshell normaliza- 
tion is enforced, using Gateaux derivatives and the set of 
Lagrange multipliers {ei}, a generalized Thomas- Fermi 
theory can be expressed in terms of partial densities and 
is equivalent to the orbital functional theory of Eqs.(8) 
[21]. 

The mathematical situation is that for more than one 
electron of either spin, the functional F s cannot be ex- 
tended to derive Euler-Lagrange equations for an un- 
structured total density p, because a Frechet functional 
derivative does not exist. However, because Gateaux 
derivatives exist, an extended functional of orbital sub- 
shell densities does exist. Independent subshell variations 
imply Euler-Lagrange equations equivalent to the OEL 
equations, valid for electronic subshell structure and for 
the exclusion principle. For noninteracting electrons, F s 
reduces to the Kohn-Sham functional T s , with the same 
mathematical properties. There is no implication that 
the density functional derivative of E xc must be a lo- 
cal potential function v xc (r), unless this is mandated by 
the simplified functional dependence assumed in the local 
density approximation (LDA). 



IV. NONLOCAL POTENTIALS 

The concept of a Gateaux functional derivative estab- 
lishes a direct relationship between orbital-indexed lo- 
cal potential functions and nonlocal potentials defined 
by orbital functional derivatives. Using exchange en- 
ergy as an example, the orbital functional derivative is 
5E x /m5(f>* — v x 4>i, and the functional differential can be 
expressed as 

5E X = ^ n i j d 3 r[S(t>*v x (f>i + cc] 



E 



, ,3 <f>*Vx<Pi 

rii I a r — — — 
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This determines a Gateaux derivative 
6E X 



Hid pi 



v xi (r), 



(13) 



(14) 



valid for arbitrary variations in the orbital Hilbert 
space. The technical problem of generalizing Kohn-Sham 
methodology to an exact theory reduces to methodol- 
ogy for such indexed local potentials. The occurrence 
of different potentials for different orbital indices can be 
compensated during self-consistent iterations by using 
off-diagonal Lagrange multipliers for orthogonalization. 
These parameters must vanish on convergence for closed- 
shell OEL equations. In the case of exchange in atoms, 
this is a variant of well-established methodology for it- 
erative solution of the Hartree-Fock equations [22]. In- 
clusion of an indexed correlation potential, derived from 
any parametrized functional E c , 



v C i(r) 



4>*Vc<t>i 



(15) 



is a straightforward generalization. 

Hartree-Fock methodology based on expansion in 
atomic Gaussian basis orbitals is valid in principle, but 
becomes impractical for large molecules and solids. For 
solids, multiple scattering theory (MST) [23] has been 
the method of choice for DFT calculations. It provides a 
common methodology for molecules, based on variational 
treatment of space-filling atomic cells [24]. In MST, the 
principal calculation over a range of orbital energies is to 
solve the OEL equations within each local atomic cell, 
or within the enclosing sphere of each such cell. Inde- 
pendent representations of a wave function within each 
cell are matched across intercell boundaries using varia- 
tional formalism [24] . Generalization of the intracell cal- 
culations to nonlocal potentials is straightforward, simply 
including an indexed parametrized exchange-correlation 
potential while integrating Schrodinger or semirelativis- 
tic equations at specified energy. The practical technical 
problem is thus reduced to the treatment of the long- 
range tails of indexed local potentials defined within each 
cell. In MST, the formalism is greatly simplified by the 
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fact that atomic basis functions are literally truncated 
at cell boundaries [24]. This means that there are no 
intercell overlap contributions. 

An approximation that may turn out to be both effi- 
cient and accurate is suggested by the MST formalism. 
This is a "local nonlocal" (LNL) model. True indexed po- 
tentials are to be used only within atomic cells, but are re- 
placed by their asymptotic forms in external cells. In gen- 
eral, these asymptotic forms arise from electrostatic mul- 
tipole potentials, and can be represented as long-range 
local potentials. If this is valid, the indexed potentials 
in each cell are augmented by an external field poten- 
tial that is summed over all other cells. As in the well- 
known Ewald expansion, this implies substantial cancel- 
lation of fields, including electrostatic screening, and may 
justify neglecting such terms during self-consistent iter- 
ations, including them as first-order perturbations after 
convergence. 

An important step toward practical methodology for 
long-range correlation has been proposed, and tested in 
an electron scattering model [25]. In a general context, 
the method is to compute, for each atomic cell, the first- 
order perturbed orbital functions due to external multi- 
pole fields. Intercell correlation energy then can be for- 
mulated in terms of these polarization pseudostates, and 
response potentials can be constructed that are correct 
at long range, but vanish smoothly at the cell origin. 
An example of dipolar response has been carried out in 
detail. The implied polarization potential is used to cal- 
culate polarization response effects in e — He scattering 
cross sections [25]. 
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